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Higher moments of net-proton multiplicity distributions are applied to 
search for the QCD critical point. In this paper, we will present mea- 
surements for kurtosis (k), skewness (S) and variance (a^) of net-proton 
multiplicity distributions at the mid-rapidity (\y\ < 0.5) and transverse 
momentum range 0.4 < px < 0.8 GeV/c for Au-|-Au collisions at ^s„„ 
= 7.7, 11.5, 39, 62.4 and 200 GeV, Cu-I-Cu collisions at = 22.4, 

62.4 and 200 GeV, d+Au coUisions at y/s^ = 200 GeV and p+p colli- 
sions at Y^s^ = 62.4 and 200 GeV. The moment products kct^ and Sa of 
net-proton distributions, which are related to volume independent baryon 
number susceptibility ratios, are consistent with Lattice QCD and Hadron 
Resonance Gas (HRG) model calculations at high energies (^s„„ = 62.4 
and 200 GeV). Deviations of kct^ and Sa for the Au-|-Au collisions at low 
energies (^s„„ = 7.7, 11.5 and 39 GeV) from HRG model calculations are 
also observed. 

PACS numbers: 25.75.Ld, 25.75.Dw 

1. Introduction 

The main goal of Beam Energy Scan (BES) program [T| at the Relativis- 
tic Heavy Ion Collider (RHIC) is to study the phase structure |2j, such as 
map the QCD phase boundary and search for the QCD critical point [3], of 
the QCD matter created in heavy ion cohision. By tuning the colliding ener- 
gies of two nuclei from ^/s^ = 200 GeV to ^s^j^ = 7.7 GeV, we can access 
various region of the QCD phase diagram. Higher moments (variance (cr^), 
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skewness (S), kurtosis (k) etc.) of conserved quantities, such as net-baryon, 
net-charge and net-strangeness, multiphcity distributions are very sensitive 
to the correlation length [H [5] and can be directly connected to the corre- 
sponding thermodynamic susceptibilities in Lattice QCD [6l |7] and Hadron 
Resonance Gas (HRG) model [8j. As the volume of the system is hard to de- 
termine, the susceptibility ratio, such as Xb^ /Xb^ and Xb^ /Xb^ i are used to 

compare with the experimental data as kcj^ = XbVxb^ ™d = XbVxb^- 
Theoretical calculations demonstrate that the experimental measurable net- 
proton (proton number minus anti-proton number) number fluctuations can 
effectively reflect the fluctuations of the net-baryon number [9]. Higher mo- 
ments analysis opens a completely new domain and provides quantitative 
method for probing the bulk properties of the hot dense nuclear matter [2]. 



2. Observables 

Experimentally, we measure net-proton number event-by-event wise, 
Np-p = Np — Np, which is proton number minus antiproton number. In the 
following, we use N to represent the net-proton number Np-p in one event. 
The average value over whole event ensemble is denoted hy fi =< N >, 
where the single angle brackets are used to indicate ensemble average of an 
event-by-event distributions. The deviation of N from its mean value are 
defined by 

6N = N- < N >= N - n. (1) 
The r^^ order central moments are defined as 

/ir-=<(<)Wr >,/Xl=0. (2) 

Then, we can define various order cumulants of event-by-event distributions 
as 



Ci = /U, C2 = H2, C3 = /is, (3) 

C„(n > 3) = - ^ f ^ _ ^ J CmfJ-n-m- (4) 
m=2 ^ ^ 

Once we have the definition of cumulants, various moments can be denoted 
as 

Then, the moments product kct^ and Sa can be expressed in term of cumu- 
lant ratio 

o C4 „ C3 
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3. Background Effects 

Therminator model [10] was applied to study the resonance decay ef- 
fect, which is a background effect for higher moment analysis, and also 
to check whether the net-proton fluctuations can reflect the net-baryon 
number fluctuations. In Fig. [H we show the event-by-event number dis- 
tributions of Au+Au 0-5% most central collisions at ^/s^ = 200 GeV 
from Therminator calculations for four cases. Fig. [T] demonstrates that 
the distribution for the net-proton with resonance decay is wider than 
the net-proton distributions without decay. By excluding the A decay nu- 
cleon, the net- (proton -|-neutron-|- A) has narrower distribution than the net- 
(proton-l-neutron) distributions. Fig. [2] shows centrality dependence of ko"^ 
of number distributions. The results for the four cases are consistent with 
each other within errors, which indicate the effects of resonance decay are 
small and the net-proton fluctuations can reflect the net-baryon fluctuations. 
The statistical errors are evaluated by Delta theorem method 



10" 

10" 



I Therminator 
Au+Au 200 GeV 
0-5% 



Net-p 
Net-(p+n) 
Net-(p+n+A) 
W/0 Decay Net-p 

0.4<p^<0.8 (GeV/c) 
|y|<0.5 




10 20 

Number 



40 



Fig. 1. Event- by-event number distributions 
for net-proton, net- (proton-t- neutron), net- 
(proton-f-neutron+A) and net-proton without 
resonance decay for 0-5% most central Au-(-Au 
collisions at ^Js^^ — 200 GeV from Therminator 
model. 
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Fig. 2. Centrality dependence of ku"^ of number 
distributions of four cases for Au+Au collisions 
at y/s^ =200 GeV from Therminator model. 



4. Experimental Method 

The data presented in this proceedings are obtained using the Solenoidal 
Tracker at RHIC (STAR). The main subsystem used in this analysis is a 
large, uniform acceptance cylindrical Time Projection Chamber (TPC) cov- 
ering a pseudo-rapidity range of \rj\ < 1 and full azimuthal coverage. To 
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ensure the purity and similar efficiency, the protons and anti-protons are 
identified with the ionization energy loss {dE/dx) measured by the TPC of 
STAR detector within 0.4 < pr < 0.8 GeV/c and mid-rapidity (|y| < 0.5). 
Centralities are determined by the uncorrected charged particle multiplic- 
ities {dNf-h/drj) within pseudo-rapidity |?7| < 0.5 measured by the TPC 
and the centrality bin width correction is used to eliminate volume fluctua- 
tions [12]. By comparing measured dNch/drj with the Monto Carlo Glauber 
model results, we can obtain the average number of participant (A^part) for 
each centrality. 



5. Results 

In this section, we present beam energy and system size dependence of 
various moments {M, a, S, k) and moment products (Sa and kcj^) of net- 
proton distributions. Those are Au+Au collisions at ^s^^ = 7.7, 11.5, 
39, 62.4 (year 2004) and 200 GeV (year 2004), Cu+Cu colhsions at ^/s^ 
= 22.4, 62.4, 200 GeV, d+Au at ^/s^ = 200 GeV (year 2003) and p+p 
collisions at ^/s^ = 62.4 (year 2006), 200 GeV (year 2009). The resuhs 
for Au-|-Au collisions at ^/s^ = 62.4 and 200 GeV have been published 
in the paper |13] . The errors shown in the figures are statistical error only. 
Centrality dependence of various moments of net-proton distributions for 



M Net-proton 

0.4<p^<0,8 {GeV/c) - 
|y|<:0.5 ^M-'' 

A''"'' 

■ ^Ci-'"' ^ « 


^ A 39 GeV 

□ 62.4 GeV 

♦ 200 GeV - 


. 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
s 

* -¥ - 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
i K 

**'».:.::,a.. 

r STAR Preliminary - 



100 200 300 100 200 300 
Npar. 

Fig. 3. Centrality dependence of various mo- 
ments of net-proton multiplicity distributions for 
Au+Au collisions at ^/s;;;; =39, 62.4, 200 GeV. 
The dashed lines shown in the figure are expecta- 
tion lines from CLT. 
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Fig. 4. Centrality dependence of various moments 
of net-proton distributions for Cu+Cu collisions 
at = 22.4, 62.4 and 200 GeV and p+p col- 

lisions at y/s^ = 62.4 and 200 GeV. The dashed 
lines shown in the figure are expectation lines from 
CLT. 



Au+Au, Cu+Cu and p+p collisions are shown in Fig. [3] and Fig. [H re- 
spectively. The M and a are found to be monotonically increasing with 
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increasing of Npart, while the S and k are decreasing. The various moments 
can be well described by the dashed lines shown in the figures, which are 
derived from Central Limit Theorem (CLT) by assuming the colliding sys- 
tem consists of many identical and independent emission sources \12\ I14j . 
Energy dependence of Sa and ko"^ for 0-5% most central Au+Au collisions 
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Fig. 5. Energy dependence of moment products (kct^ and Sa) of net-proton distri- 
butions for 0-5% most central Au-|-Au collisions. The red dashed lines denote the 
HRG model calculations, and the empty markers denote Lattice QCD results [5]. 



are shown in Fig. [5j We find that the data are consistent with Lattice QCD 
and HRG model calculations at high energies (^s^^ = 62.4 and 200 GeV), 
while deviations from HRG model are observed at low energies {y/s^ = 
7.7, 11.5 and 39 GeV). The possible reasons for the deviations are discussed 
in da [IS]. 

6. Summary and Outlook 

Higher moments of net-proton distributions are applied to search for 
the QCD critical point and probe the bulk properties of QCD matters. In 
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summary, we present the measurements of higher moments of net-proton 
distributions for Au+Au, Cu+Cu, d+Au and p+p coUisions from STAR 
experiment. The moment products ko"^ and Sa of net-proton distributions 
from 0-5% most central Au+Au coUisions are consistent with Lattice QCD 
and HRG model calculations at high energies {^/s^ = 62 A and 200 GeV), 
while the results are smaller than HRG model calculations at low energies 
iV^NN ~ ^ 11-5, 39 GeV). The analysis of data from another two energies 
at ^/s^ = 19.6 and 27 GeV, which were collected in the year 2011, are 
ongoing. 
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